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ABSTRACT 

Binary stars hosting exoplanets are a unique laboratory where chemical tagging can be performed to measure with 
high accuracy the elemental abundances of both stellar components, with the aim to investigate the formation of 
planets and their subsequent evolution. Here, we present a high-precision differential abundance analysis of the XO-2 
wide stellar binary based on high resolution HARPS-N@TNG spectra. Both components are very similar K-dwarfs 
and host planets. Since they formed presumably within the same molecular cloud, we expect they should possess the 
same initial elemental abundances. We investigate if the presence of planets can cause some chemical imprints in the 
stellar atmospheric abundances. We measure abundances of 25 elements for both stars with a range of condensation 
temperature Tc ~ 40 — 1741 K, achieving typical precisions of ~ 0.07 dex. The North component shows abundances 
in all elements higher by -1-0.067 ± 0.032 dex on average, with a mean difference of 4-0.078 dex for elements with 
Tc > 800 K. The significance of the XO-2N abundance difference relative to XO-2S is at the 2(j level for almost all 
elements. We discuss the possibility that this result could be interpreted as the signature of the ingestion of material 
by XO-2N or depletion in XO-2S due to locking of heavy elements by the planetary companions. We estimate a mass 
of several tens of M 0 in heavy elements. The difference in abundances between XO-2N and XO-2S shows a positive 
correlation with the condensation temperatures of the elements, with a slope of (4.7 ±0.9) x 10“® dex K“^, which could 
mean that both components have not formed terrestrial planets, but that first experienced the accretion of rocky core 
interior to the subsequent giant planets. 
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* Based on observations made with the Italian Telescopio 
Nazionale Galileo (TNG), operated on the island of La Palma 
by the INAF - Fundacion Galileo Galilei at the Roche de los 
Muchachos Observatory of the Instituto de Astroflsica de Ga- 
narias (lAC) in the framework of the large programme Global 
Architecture of Planetary Systems (GAPS; P.I. A. Sozzetti). 


1. Introduction 


A well-established direct dependence exists between 
the occurrence rate of giant planets and the metal 


content of their main-sequence hosts 

Gonzaled 1997: 

Butler et al. 20001: Laws et al.l I200S 

: Santos et al.l 

2004 iFischer & Valentil 20051: 

Udrv & SantosI 20071: 

Sozzetti et al. 20091: Mortier et al.l 

201211. 

Two main see- 


narios have been invoked to explain the observational 
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evidence. On the one hand, the enhanced likelihood of 
forming giants planet in higher-metallicity disks can be 
seen as a direct con s equence of the c o re-accretion scenari o 
(|Pollack et al.l [To^ llda fc LinI l2004 lAlibert et al.ll200^ . 
because the expected larger reservoirs of dust grains 
in more metal-ric h environments (e.g., iPei fc Falll ll^ 
iDraine et al.l 120071) will favor the build-up of cores that 
will later accrete gas. In this case, the observed giant 
planet - metallicity connection has a primordial origin. 
On the other hand, it is possible that the outcome of the 
planet formation process might directly affect the chemical 
composition of the host stars. Two different versions of 
this second scenario have been proposed. In the first 
one, disk-planet or planet-planet interactions may cause 
the infall of me tal-rich ("rocky") material onto the star 
(|Gonzaled Il99i11 . In this case, the presence of planets is 
expected to positively correlate with metallicity. A similar 
effect is thought to be in action in the Solar System 
(|Murrav et al.l l200lh . where however the total mass of 
rocky material accreted by the Sun after the shrinking 
of its convective zone 0.4 M 0 ) is possibly too low for 
causing an appreciable variation of its chemical composi¬ 
tion. Alternatively, the formation of planetary cores may 
prevent rocky planeta ry material to be accreted by the star 
l|Melendez et al.l[2009jl . In this case, the presence of planets 
should negatively correlate with metallicity, in the sense 
that depletion of refractory elements relative to volatiles 
should be observed. Some evidence that t his might be the 
case f o r the Sun has bee n presented by iMelendez et al.l 
l|2009ll . IChambersI (j201Clll estimated that the deficiency 
of refractory elements found by these authors would be 
canceled if ^ 4 AI^ of rocky material wo uld be added to 
the pr esent-day solar convective envelope. IMelendez et al.l 
(|2009^ also concluded that solar twins without close-in 
giant planets chemically resemble the Sun, suggesting that 
the presence of such planets might prevent the formation 
of Earth-size planets. Of course, it is possible that all 
these (and other) mechanisms are at work, rendering 
the discernment of thei r relative roles a rather complex 
undertaking (see, e.g., [Gonza lez H ernandez et all 120131 : 
lAdibekvan et al.ll2014l : lMaldonado et al.ll20r5h . 


While accretion of plane tary material may caus e huge 
effects in white dwarfs le.g.. IZuckerman et al.ir2003n . vari¬ 
ations of surface abundances are not expected to be 
large in FGK main sequence stars due to the dilu¬ 
tion of the material wit hin the outer convective enve¬ 
lope (jMurrav et al.l 1200 ill . In general, the most promi¬ 
nent effect should be an over- or under-abundance of ele¬ 
ments as a function of the grains condensation temperature 
(Tn). which is different for different elements ([Smith et al.l 
l2001h . Various authors have tried to show that such 
anomalies may be present in stars with p l anets (see, e.g., 
Smith et al. 2001: Laws fc Gonzalez 200lt Melen dez et al'l 


20091 : ^mirez et al.ll201(3r Maldonado et al.ll2015ll . Reveal- 

ing these trends requires the development of rather inge¬ 
nious techniq ues of very accurate differential abundance 
analysis (e.g., iRamirez et al.l [201411 . In spite of these ef¬ 
forts, the results, if any, are still quite elusive. Obtaining 
accurate differential abundance s is easier in the case of 
stars member s of clusters (e.g., Yong_et_aI, 2^^ and bi- 


nary 1 


systems (Grattpnetak 2Q0lJ: LawsfcCjon^alej 2Q01 


Desidera et al.l20M . I2006I : iRamfrez et al.l 20111 iLiu et al.l 

201411 . because here many obseryational uncertainties re¬ 
lated to distances and reddening can be considered 


common-mode effects. In addition, the most accurate re¬ 
sults are obtai ned comparing stars that are yery similar 
to each other (iGratton et al-l 120011 : iDesidera et al.l l2004 
iLiu et aDl2014ll . At the same time, it is also crucial that 
such stellar samples be accurately scrutinized for the pres¬ 
ence of planets around them. 

XO-2 is a wide yisual binary (separation 31 arc- 
sec, correspondin g to a projected distance of ~ 4600 AU; 
iBurke et al.ll20f)7TI compose d of two uneyolyed st ars of yery 
similar masses (^ 0.97M^: iDamasso et al.l[201.4) . and with 
the same proper motions (LSPM J0748-I-5013N and LSPM 
J0748+5013S in lLeoine fc Shar^l2005ll . Both stars have no 
detectable a,mounts of lithium, s upporting the old as e of 
the system (jPesidera et al.l[20l4 iDamasso et al.ll201jR ). A 
transiting g i ant pl anet was discovered around XO-2N by 
iBurke et al. l|2007ll . with a period P = 2.616 d, a mass 
Mp = 0.597 Mj, and a semi-major axis a = 0.037 AU. Re¬ 
cently, XO-2S was also shown to host a planetary system 
composed of at least two planets (Mp sinz = 0.259 Mj, 
P = 18.157 d, a = 0.134 AU and M„ sini = 1.37 0 Mj, 
P = 120.80 d, a = 0.476 AU; IDesidera et al.l[2?)l4 here¬ 
after Paper I). In a subsequent characterization work on 
the XO-2 system, the existence was established of a long¬ 
term curvature in the radial velocity (RV) measurements of 


XO-2N, suggesting the possible prese nce of a second planet 
in a much wider orbit with P > 17 yr ([Damasso et al.l[2015L 
hereafter Paper II). The reason for the difference in the ar¬ 
chitectures of the planetary systems around the two stars 
is presently unclear. 

An accurate differential analysis of the iron content 
in the two components of XO-2 revealed a small but 
clearly significant difference of 0.054 ± 0.013 dex, XO-2N 
being more iron-rich than XO-2S (Paper II). This is a 
rare occurrence among binary systems, as in most cases 
the two components have v ery similar metal abundances 
(jPesidera et al.l l2004 l200(^. This resu lt see ms to be in 
contrast with the work by iTeske et al.l l|2013ll . who found 
the two stars similar in physical properties and in iron 
and nickel abundances, with some evidence of enrichment 
in C and O abundance o f XO-2 N compared to XO-2S. 
Very recently, iTeske et al.l l|2015ll explored the effects of 
changing the relative stellar temperature and gravity on 
the elemental abundances. Using three sets of parameters, 
they found that some refractory elements (Fe, Si, and Ni) 
are most probably enhanced in XO-2N regardless of the 
chosen parameter^ Therefore, in the XO-2 system, sim¬ 
ilar to other binary systems hosting planets, the results 
are still ambiguous. In fact, there is stil l no consensus on 
the d i fferences between 16CygA and B (iLaws fc Gonzal^ 
20W; Takedal I 2 OO 5 I: fSchuler et al.l 2011a: Ramirez et al 


2011 ; Metcalfe et al.l 120121 : iTucci Maia et al.l 12014 , and 
references therein), while no significant abundance dif- 


^ In iDamasso et aP (1201511 . we measured the XO-2 rotation pe¬ 
riods, finding a value of 41.6 days in the case of XO-2N and 
a value between 26 and 34.5 days for XO-2S, combined with 
slight decreases in the activity levels of both stars at the end 
of the 2013-2014 season. Analyzing new APACHE photometric 
data taken from October 2014 to April 2015, both these findings 
are not strengthened, probably as a result of the still decreased 
magnetic activity. 

^ After the present pap er was submitted, an independent work 
by IRamfrez et all (120151 1 based on different spectroscopic data 
and methodologies was published, with good agreement and con¬ 
sistent results with our analysis. 
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ferences be tween componen ts in the singl e-host systems 


HAT -P-l (iLiu et al.l ^oTl. HD 13256 3 (iDesidera et all 


1201 Ih . HD 106515 I Desidera et al.l 201^ and the dual-host 
HD 20782/1 system ( Mack et al.ll20lj^ were found. 


The XO-2 binary system is of particular interest. Be¬ 
ing composed of two bright, near-twin stars with well- 
characterized planetary systems around them, it is one of 
only four kno wn dual-planet-hos ting binaries, t ogether with 
HD 2 0782/1 iTones et al.l l200 6ll. Kenler 132 (iR.owe et all 
[20 i 3), and WASP 94 I Neveu-Vanmalle et al.l[^014^ . As we 
have mentioned above, planet-hosting binaries provide use¬ 
ful laboratories to try to decouple the effects of the stel¬ 
lar system initial abundances due to the star forming cloud 
from the possible effects of the planetary system on the star 
chemical composition. It is therefore important to perform 
a more extensive abundance analysis of the XO-2 system, 
not only to confirm the abundance difference we found for 
iron between two components, but also to establish if there 
is any correlation of the abundance differences with the 
condensation temperatures or other quantities. 


As a follow-on of Papers 1 and 11 and within the context 
of the GAP£ ^ (Gl obal Architecture of Planetary Systems; 


ICovino et al 


201311 programme, we present an accurate dif¬ 


ferential abundance analysis of 25 elements for the two com¬ 
ponents of the XO-2 system (XO-2N and XO-2S). The use 
of the procedure described in Paper 11 allows us to derive 
the differences in elemental abundances (and stellar param¬ 
eters) with very high accuracy and unveil small differences 
among components, removing several sources of systematic 


The outline of this paper is as follows. We first briefly 
present in Sect. [3] the spectroscopic dataset. In Sect. [31 we 
describe the differential analysis of the elemental abun¬ 
dances for 25 species. We then discuss the behavior of 
the elemental abundance differences with the condensation 
temperature (Sect.S]) and the origin of these differences 
(Sect.|S|). In Sect. ini our conclusions are presented. 


2. Observations 

We observed both XO-2 components with the high resolu- 
tion HARPS-N@TNG {R ~ 115 000, A ~ 3 900 - 6 900 A; 
ICosentino et al.ll2012h spectrograph between November 20, 
2012 and October 4, 2014. Solar spectra were also obtained 
through observations of the asteroid Vesta. The spectra re¬ 
duction was obtained using the 2013 November version of 
the HARPS-N instrument data reduction software (DRS) 
pipeline. A detailed description of the observations and data 
reduction is reported in Paper H. 

As done in Paper H, we measured the elemental abun¬ 
dances co-adding the spectra of both components and 
Vesta, properly shifted by the corresponding radial veloc¬ 
ity, to produce merged spectra with signal-to-noise ratios 
(S/N) around 300 (for the targets) and 400 (for the aster¬ 
oid) per pixel at A ~6000 A. 


® HD 20781 should host two Neptune-mass planets, as re¬ 
ported by Mayor et al. (2011), but the discovery was not con¬ 
firmed /published. 

* http://www.oact.inaf.it/exoit/EXO- 
IT/Projects/Entries/2011/12/27_G APS.html 


3. Spectral analysis 

3.1. Differential elemental abundances 

Elemental abundance^ were measured using the differen¬ 
tial spectral analysis described in Paper H, which is based 
on th e pr escriptions wid e ly describe d in iGratton et al.l 
(|200lfl and IDesidera et ^ l|2004 I^OOfl . In particular, this 
method gives very accurate results when the binary compo¬ 
nents are very similar in stellar parameters, as in the case of 
the XO-2 system. We thus considered an analysis based on 
equivalent widths (EWs), where a strict line-by-line method 
was applied, for both stars and the solar spectrum. Then, 
XO-2S was analyzed differentially with respect to the N 
component, which was used as a reference star. This al¬ 
lowed us to obtain the differences in stellar parameters and 
elemental abundances between XO-2N and XO-2S. This ap¬ 
proach minimizes errors due to uncertainties in measure¬ 
ments of EWs, atomic parameters, solar parameters, and 
model atmospheres. Moreover, being the two components 
of this binary system very similar to each other, this reduc es 
the possibility of systematic errors (jPesidera et al.|[2006ll . 

We considered as stellar parameters (effective temper¬ 
ature Teff, surface gravity logt/, microturbulence ^) and 
iron abundances ([Fe/H]), together with t heir uncertain- 
ties, those reported in our Paper H, where iKurucd (|l993ll 
grids of plane-paralle l model atmospher es were used. We 
refer to that work and iBiazzo et al.l (l2012fl for a detailed de¬ 
scription of the method to measure abundances (and their 
uncertainties). Once stellar parameters and iron abundance 
were set, we computed elemental abundances of several re¬ 
fractories, siderophiles, silicates, and volatiles elements (ac- 
cording to their cosmoche mical character; see, e.g.. Fig. 4 of 
Anders fc Grevess'3 1 198911 using the MOOG code (ISiiedenI 
19731 version 20141 and the abjind and blends drivers for the 
treatment of the lines without and with hyperfine structure 
(HFS), respectively. We therefore obtained the abundance 
of 25 elements: G, N, O, Na, Mg, Al, Si, S, Ga, Sc, Ti, V, 
Gr, Mn, Fe, Go, Ni, Cu, Zn, Y, Zr, Ba, La, Nd, and Eu. 
For three elements (Ti, Gr, and Fe) we measured two ion¬ 
ization states, while for the other elements only one species 
(first or second ionization state) was measured. As done in 
Paper H for the iron lines, here the EWs of each spectral 
line were measured on the one-dimensional spectra interac¬ 
tively using the splot task in IRAE0. The location of the 
local continuum was carefully selected tracing as much as 
possible the same position for each spectral line of both bi¬ 
nary system components and the asteroid Vesta. This was 
done with the aim to minimize the error in the selection 
of the continuum. We also excluded features affected by 
telluric absorption. 

We compi led a long list of spe ctral lines, starting from 
the work by iBiazzo et al.l (I2fll2fl . In particular, our line 
list comprises Nai, Mg I, All, Sii, Cal, Til, Tin, Cri, 
Cm, Fei, Fell, Nil, and Zni, complemented with ad¬ 
ditional lines and atomic parameters for Nai, All, Sii, 
Til, T ill, Cri, Nil, and Zni taken from ISchuler et al.l 
(l20m and lSozzetti et ahl (120061) to increase the statistics. 

® Throughout the paper, the abundance of the X element is 
given in log -|- 12. 

® IRAF is distributed by the National Optical Astronomy Ob¬ 
servatory, which is operated by the Association of the Universi¬ 
ties for Research in Astronomy, Inc. (AURA) under cooperative 
agreement with the National Science Foundation. 
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Wavelength (A) Wavelength (A) 


Fig. 1. Comparison between observed (black line) and syn¬ 
thetic spectra for XO-2S in the region around the [O l] line 
at 6300.3 A. The synthetic spectra have been computed for 
logn(O) = 9.15 (blue line), 9.20 (red line), and 9.25 (green line), 
respectively. Note that several CN lines are also present in this 
spectral region; due to coupling between C and O, the strength 
of these features is anti-correlated with that of O. 


In seven cases (Cl, Si, Sen, Vi, Mni, Coi , and Cul), we 


conside re d the line lists by Kuni^ (119931). Schuler et al.l 


(|2011ljl . I.Tohnson et al. ( 2006117 and IScott e t al. (l2015bD . 


where the HFS bv iJohnson et al.l (I2006D and Kurued (1199^ 
wa s adopted for Sc, V, Mn, Cu , and Co. Solar isotopic ratios 
by lAnders fc Crevessg (|l989h were considered for Cu (i.e. 
69.17% for 63Cu and 30.83% for ^^Cu). For the s-process 
elements Yii, Zrii (first peak) and Ban, Laii (second 
peak) and the r-proces s elements Ndi i and Full we con¬ 
sidere d the line lists by I Johnson et al.l (12006 ] , iLiung et al.l 


(1200611. i Prochaska et a.i.1 (|200^)^ . iLawler et ^ (|2001h ), and 
iDen Hartog et a l.l (120031) where t he HFS by lGallagher et all 
(|2010ll and iLawler et al.l (|2001all was adopted for Ba and 
La, re spectively. Solar isotopic ratios bv lAnders fc Crevessel 
l|l989ll were considered for Ba (i.e. 2.417% for 7.854% 

for i36Ba, 71.70% for ^^SBa, and 6.592% for ^^SBa, 11.23% 
for i37Ba) and Eu (i.e. 47.8% for i^iEu and 52.2% for 

i53Eu). 


CNO abundances were measured as follows. First, O was 
measured by u sing the synth driver within the MOOG code 
(|Snedenll 19731 version 2014) fitting synthetic spectra to the 
zone around the [Oi] line at 6300.3 A (see Fig. [T]). Then 
adopting the estimated O abundance, C was measured from 
features due to the C 2 Swan system at 5150 ((0,0) band) 
and 5635 A ((0,1) band: see Fig. [5]). Finally, N was measured 
from several CN features at ^4195 and ^4215 A, adopting 
the measured C and O abundances. The process was re¬ 
peated for both the stars and for the Vesta spectrum until 
convergence of the derived values. Note that the procedure 
was strictly differential, using the same features in each of 
the stars. Line lists were built using the latest version of 


Fig. 2. Comparison between observed (black line) and syn¬ 
thetic spectra for XO-2S in the region around the bandhead 
of C 2 at 5635 A. The synthetic spectra have been computed for 
logn(C) = 8.83 (blue line), 8.88 (red line), and 8.93 (green line), 
respectively. 


seQ _ 

lists for CN are f rom lSneden et ah I (|20T1 a nd C 2 are from 


the VALD databas qll for the atomic lin es, while the line 


iR.am et ^ l|2014^ and Brooke et al. ( 201311 . Carbon abun¬ 
dances were also derived from atomic C i. Since C abun¬ 
dance was derived with both EWs and spectral synthesis, 
as final abundance we considered the mean value coming 
from these two methods. 

Table [T] lists the solar elemental abundances we ob¬ 
tained using spectra of the Vesta asteroid and the follow¬ 
ing solar parameters: Tpff c^ = 5770 K, logq,j^ = 4.44, and 
^0 = 1.10 km/s (see, e.g.. lBiazzo et al]l2012ll . In the same 
table, the most recent determinations of solar abundances 


bv lAsplund et al.l (120091) . Scott et ^ l|2015al) . IScott et ^ 
(|2015bl) . and Crevesse et al. ( 20151) are also given. We stress 
that the latter values were obtained using 3D+NLTE mod¬ 
els. Our determinations are therefore in good agreement 
with these recent literature values. 


In Table [5] we present the final elemental abundances 
(relative to the solar ones) and the differential abundances 
(N—S) of all elements, together with their respective er¬ 
rors. The estimate of the uncertainties on each elemental 
abundance and on differential [X/H] was done as widely 
described in Paper II (see also references therein). Table |3] 
lists the mean internal errors in abundances determination 
due to uncertainties in stellar parameters for both com¬ 
ponents, where the stellar parameters and iron abundance 
are those obtained and reported in Paper II. In the end, the 
differences in parameters (in the sense of N minus S compo¬ 
nent) are: ATgff = —35±8K, Alogg = -|-0.010±0.020 dex, 
A^ = -0.07 ± 0.07 km s"!, and [Fe/H]= -b0.055 ± 0.013 
dex. Our analysis demonstrates that: i) the S component is 
slightly hotter than the N component, as reported in Paper 


^ http://vald.astro.uu.se/ 
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II; a) besides the enhancement in iron found in our Pa¬ 
per II, the XO-2 binary system is more rich than the Sun 
in almost all elements (see Fig. |31 where the differential 
abundances of the N minus S component are displayed as 
a function of the atomic number Z and the ionization po¬ 
tential /); in) XO-2N is enhanced in all refractory elements 
and less enhanced in average in volatiles (see Fig. 01 where 
the differential abundances are displayed as a function of 
the condensation temperatures. Hin ts of a similar behavior 
were reported bv lTeske et al.l (|2015!l for Fe, Si, and Ni. 

All differential elemental abundances have errors in the 
range ^ 0.01 — 0.05 dex (with the exception of the eu¬ 
ropium), which further demonstrates the advantages of a 
strictly differential analysis. The mean abundance differ¬ 
ence of all elements between XO-2N and XO-2S is A[X/H]= 
+0.067 ±0.032 dex, with a mean difference of 0.025 dex for 
the volatile elements (T^ < BOOK) and of ^ 0.078dex for 
the elements with Tq > 800 K. For the neutron-capture el¬ 
ements (e.g., Yii, Zrii, Ban, Laii, Ndii, Full) the mean 
difference is 0.087 dex. No elemental abundance of XO-2N 
differs from that of XO-2S by more than ^0.11 dex, and all 
but three species (namely, N, S, and Eu) have differences 
detected at ^ 2cr level. 


Table 1. Gomparison between our measured solar abundances 
and standard values from the literature. 


Species log riThis work It^S ^Literature Reference 


c 

8.50 

± 

0.03 

8.43 

± 

0.05 


N 

7.93 

± 

0.04 

7.83 

± 

0.05 

(1) 

0 

8.93 

± 

0.05 

8.69 

± 

0.05 

(1) 

Na 

6.305 

± 

0.007 

6.21 

± 

0.04 

(2) 

Mg 

7.590 

± 

0.044 

7.59 

± 

0.04 

(2) 

Al 

6.490 

± 

0.115 

6.43 

± 

0.04 

(2) 

Si 

7.558 

± 

0.052 

7.51 

± 

0.03 

(2) 

S 

7.210 

± 

0.014 

7.12 

± 

0.03 

(2) 

Ga 

6.347 

± 

0.035 

6.32 

± 

0.03 

(2) 

Sc 

3.135 

± 

0.049 

3.16 

± 

0.04 

(2) 

Til 

4.960 

± 

0.043 

4.93 

± 

0.04 

(3) 

Ti II 

4.982 

± 

0.048 




V 

3.856 

± 

0.039 

3.89 

± 

0.08 

(3) 

Gri 

5.662 

± 

0.033 

5.62 

± 

0.04 

(3) 

Gr II 

5.655 

± 

0.049 




Mn 

5.407 

± 

0.064 

5.42 

± 

0.04 

(3) 

Fei 

7.530 

± 

0.049 

7.47 

± 

0.04 

(3) 

Fell 

7.535 

± 

0.053 




Go 

4.894 

± 

0.074 

4.93 

± 

0.05 

(3) 

Ni 

6.279 

± 

0.048 

6.20 

± 

0.04 

(3) 

Gu 

4.213 

± 

0.074 

4.18 

± 

0.05 

(4) 

Zn 

4.553 

± 

0.025 

4.56 

± 

0.05 

(4) 

Y 

2.165 

± 

0.051 

2.21 

± 

0.05 

(4) 

Zr 

2.600 

± 

0.031 

2.59 

± 

0.04 

(4) 

Ba 

2.235 

± 

0.007 

2.25 

± 

0.07 

(4) 

La 

1.105 

± 

0.007 

1.11 

± 

0.04 

(4) 

Nd 

1.460 

± 

0.044 

1.42 

± 

0.04 

(4) 

Eu 

0.543 

± 

0.042 

0.52 

± 

0.04 

(4) 


No tes: (1) Asplu nd et al.l (l2009f): 12) Scottet_alJ l|2015bll : 
('Si lScott et al. 120 15 all : 141 iGrevesse et al. 1 201511 . 


Table 2. Differential elemental abundances for the XO-2 sys¬ 
tem. 


Element 

rjia 

XO-2N'’ 

XO-2S'’ 

A[X/H]" 

[C/HJ 

40 

0.41 

± 

0.04 

0.36 

± 

0.05 

0.050 ± 0.02 

[N/H] 

123 

0.47 

± 

0.05 

0.49 

± 

0.05 

-0.020 ± 0.02 

[O/H] 

180 

0.32 

± 

0.03 

0.27 

± 

0.03 

0.050 ± 0.02 

[Na/H] 

958 

0.485 

± 

0.043 

0.445 

± 

0.043 

0.040 ±0.011 

[Mg/H] 

1336 

0.377 

± 

0.051 

0.333 

± 

0.048 

0.044 ±0.011 

[Al/H] 

1653 

0.440 

± 

0.139 

0.343 

± 

0.141 

0.097 ±0.017 

[Si/H] 

1310 

0.471 

± 

0.082 

0.419 

± 

0.082 

0.052 ±0.010 

[S/H] 

664 

0.530 

± 

0.020 

0.525 

± 

0.025 

0.005 ± 0.027 

[Ca/H] 

1517 

0.365 

± 

0.075 

0.301 

± 

0.065 

0.064 ± 0.023 

[Sc/H] 

1659 

0.462 

± 

0.067 

0.388 

± 

0.060 

0.074 ± 0.009 

[Ti/H]* 

1582 

0.456 

± 

0.096 

0.351 

± 

0.091 

0.105 ±0.030 

[V/H] 

1429 

0.431 

± 

0.060 

0.329 

± 

0.062 

0.102 ±0.013 

[Cr/Hj* 

1296 

0.423 

± 

0.059 

0.339 

± 

0.054 

0.084 ± 0.043 

[Mn/H] 

1158 

0.576 

± 

0.065 

0.506 

± 

0.074 

0.070 ± 0.043 

[Fe/H]* 

1334 

0.370 

± 

0.072 

0.315 

± 

0.078 

0.055 ±0.013 

[Co/H] 

1352 

0.524 

± 

0.127 

0.448 

± 

0.126 

0.076 ± 0.009 

[Ni/H] 

1353 

0.538 

± 

0.094 

0.455 

± 

0.096 

0.083 ± 0.020 

[Cu/H] 

1037 

0.627 

± 

0.109 

0.542 

± 

0.107 

0.085 ±0.019 

[Zn/H] 

726 

0.467 

± 

0.044 

0.427 

± 

0.059 

0.040 ± 0.023 

[Y/H] 

1659 

0.472 

± 

0.134 

0.360 

± 

0.110 

0.112 ±0.036 

[Zr/H] 

1741 

0.460 

± 

0.084 

0.350 

± 

0.084 

0.110 ±0.026 

[Ba/H] 

1455 

0.410 

± 

0.022 

0.310 

± 

0.049 

0.100 ± 0.049 

[La/H] 

1578 

0.530 

± 

0.036 

0.460 

± 

0.022 

0.070 ± 0.014 

[Nd/H] 

1602 

0.547 

± 

0.151 

0.460 

± 

0.147 

0.087 ±0.013 

[Eu/h] 

1356 

0.520 

± 

0.108 

0.477 

± 

0.104 

0.043 ± 0.091 


Notes: “ As condensation temperatures w e considered the 
50% Tc values derived by iLoddersI (l2003ll at a total pres¬ 
sure of 10“'^ bar. ^ Abundances relative to the Sun. Abun¬ 
dances of XO-2N relative to XO-2S. Errors refer to the dif¬ 
ferential elemental abundances (see text). * Ti, Cr, and Fe 
abundances refer to the values obtained with the first ion¬ 
ization state. 


3.2. Stellar population of the XO-2 components 

Elemental abundances, together with Galactic kinematics, 
are often used to identify a star as belonging to a spe¬ 
cific stellar population. This is because the thin disk and 
the thick disk within the Milky Way are discrete popula¬ 
tions showing distinct chemical distributions in several ele¬ 
ments, and in particular in a-elements. In fact, the Galactic 
thin and thick disks appear to overlap significantly around 
—0.7 ^ [Fe/H] ^ 0.0 dex when the iron abundance is used 
as a reference, while they are separated in [a/Fe] (see, e.g., 
iBensbv et al.ll200^ iMishenina et al.l 120041) . 

Here, we investigate to which population the XO- 
2 stellar system belongs to taking advantage o f 431 
stars within the catalogue by [Soubiran|^^Bira2_|200^ 
and 9 06 FGK dwarf stars analyzed by lAdibekvan et al.l 
(|2012fl with thin-disk, thick-disk, and halo member¬ 
ship probabilities (Pthin, Pthick, Phaio) higher than 90%. 
We use as abundance of a-elements the definition 
[a/Fe] = |;([Mg/Fe]+[Si/Fe]+[Ga/Fe] + [Ti/Fe]), since mag¬ 
nesium, silicon, calcium, and titanium are the common el¬ 
ements analyzed by those authors and by us, and also be¬ 
cause they show similar behavior in the chemical evolution 
within our Galaxy. Figure[S]displays the position of the XO- 
2N and XO-2S stars in the [a/Fe] versus [Fe/H] diagram. 
Based on this sort of “chemical indicator”, the stellar system 
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Fig. 3. Elemental abundance difference between XO-2N and 
XO-2S as a function of the atomic number Z {upper panel) and 
ionization potential I {lower panel). In both panels, dotted line 
refers to null difference, dashed line represents the mean value 
of the difference, and the filled grey area is the ±1 ct level. 



Fig. 5. [a/Fe] versus [Fe/H] for the lAdibekvan et all (I2012I I 
sample, where open circles, open squares, and asterisks repre¬ 
sent stars with Fthin > 0.90, Pthick > 0.90, and Phaio > 0.90 
ffrom iBensbv et aLl200,'lf), respectively. Filled circle s and squares 
show the stars within the lSoubiran fc G irardI ll 2005fl sample with 
Pthin > 0.90 and Pthick > 0.90 ffrom lMishenina et al.ll200J'l . re¬ 
spectively. The positions of the two targets XO-2N and XO-2S 
are displayed with big yellow and green five-pointed star sym¬ 
bols, respectively. 



Tc (K) 


Fig. 4. Differential elemental abundances of XO-2N — XO-2S 
versus condensation temperatures. Solid blue line is the un¬ 
weighted linear least-square fit to our data. The horizontal dot¬ 
ted blue line represents the average of the volatiles, while the 
other dotted blue line is the trend of the refracto ry elements. The 
dashe d red line is the mean trend obtained bv (Melendez et aP 
ll2009ll for eleven solar twins with respect to the Sun, after a 
vertical shift was applied to match the refractory elements. 


studied in this work appears to belong to the Galactic thin 
disk population. This resul t gives support to previous find¬ 
ings based on kinematics (iBurke et al.ll200^ and analysis 
of the Galactic orbits (Paper II). 


4. Abundance differences versus condensation 
temperatures 


Possible evidences of planetesimal accretion can be inves¬ 
tigated by searching for any dependence between elemen¬ 
tal abundances and condensation temperatures. This is be¬ 
cause any accretion event would occur very close to the 
star (i.e. in a high-temperature environment). Since refrac¬ 
tory elements (e.g., a- and Fe-peak elements) condense at 
high Tc, they might be added in large quantities with re- 
spe ct to volatiles eleme nts (e.g., C, N, O, Zn) with lower 
Tc ((Sozzetti et al.ll^OOfil and references therein). 

Any trend of [X/H] with Tc can be quantified in 
terms of a positive slope in a linear fit. The pos¬ 
sibili ty of a trend of [X/H] with Tc llSozzetti et al.l 


20061: Gonzalez Hernandez et al.l (20131 : (Liu et al.l (2014 
Mack et al.[ l20lf and references therein) and its depen¬ 


dence on stell ar age, surface gravity, and mean galactocen- 
tric distance ((Adibekvan et al][2014r) has been investigated 
by several authors, and often no statistically convincing 
trend could be found. 


The similarity in stellar parameters between the XO-2 
binary components allowed us to derive differential elemen¬ 
tal abundances with a high level of confidence. Figure [4] 
shows the differential abundances of the XO-2N and XO- 
2S stars relative to each other versus Tc in the range 40- 
1741 K (see also Table [5]). We remark the importance to 
have as much as possible an extended range of condensa¬ 
tion temperatures to put constraints in such a trend. The 
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slope of the unweighted linear least-square fit is positive, 
with a value of (4.7 ± 0.9) x 10“® dex K“^. We calcu¬ 
lated the Pearson’s r correlation coefficient by me ans of 
the IDL0 procedure CORRELATE (iPress et al.lll986ll . This 
test indicates that A[X/H] and Tc are positively corre¬ 
lated, with r ~ 0.72. Considering only the elements with 
higher values of Tq (i.e. > 800 K), the slope and r become 
(6.1 ±2.1) X 10“^ dex and 0. 57, respectively. More¬ 
over, using Monte Carlo simulations (jParke Lovd fc Francd 
l2014^ , we found that the probability that uiicorrelated ran¬ 
dom datasets can reproduce the observed arrangement of 
point is only ~ 0.9%. This proves that the observed positive 
correlation between the differential elemental abundances 
of XO-2N related to XO-2S and the condensation tempera¬ 
tures is significant. A similar peculiar positive trend in the 
observed differential abundances with Tq has been found 
in solar twin stars known to have close-in giant planets 
(see Fig. il), while the majority of solar analogues without 
revealed giant planets in radial vel ocity monitoring show 
a negative trend, such as the Sun ([Melendez et al.l 12009(1 . 
This result would seem to be consistent with the expecta¬ 
tion that the presence of close-in giant planets p revents the 
forma tion and/or survival of terrestrial planets (llda fc LinI 

I2OOI . 


4.1. Comparison with previous works 

In Paper II we have found that XO-2N is richer in Fe 
abundance than the S component. Here, we confirm this 
finding for many more elements (in particular, the refrac¬ 
tories). This result seems to be consistent with the re¬ 
cent work by iTeske et al.l ([201^, who found that the N 
component is most probably enhanced in Si, Fe, and Ni 
when compared to XO-2S. In the case of I6Cyg A/B, 
no consensus was reached, as some auth ors find differ¬ 
ences in the differential abu ndances fe.g.. I^amfrez et al.l 
l2011f iTucci Mala et all 1201411 . and others claim the stars 
are chemically homoge neous (with the exception of Li and 
B). iMack et al.l (|201# find non-zero differences for O, Al , 
V, Cr, Fe, and Co of HD20781/2. Again. [Liu et al.l l|2014ll 
found A[X/H] yalues consi stent with zero in the HAT-P- 
1 G0±F8 binary. Recently, iTeske et al.l (l2015f) claim that 
perhaps we are seeing an effect in A[X/H] trends due to 
Teff, moving from hotter stars with no differences to cooler 
stars with more significant differences. They also warn that 
we have to be cautious about this question, as it is hard at 
the present day to have an answer due to the different anal¬ 
ysis methods and the host star nature of these systems (i.e. 
wheth er or not they host p lanets). Moreover, as pointed 
out bv ISozzetti et al.l (l2fl06li . there is an intrinsic difficulty 
in determining with accuracy elemental abundances for a 
large set of elements in cool stars without the danger to 
introduce greater uncertainties in the results. 

Simila rly to the analysis of several binary systems 
(see, e.g.j Rmnfrgzjd^jl 200 ^ 2^9 IT Schulgrjdjl 201^ 


iLiu et all I20l4 iMack et al.f 2014 ; Tucci Maia et al.ll201^ 

we find a positive trend with Tq for the elemental abun¬ 
dances. Our results are consistent with accretion events of 
H-depleted rocky material on both components. In single 
stars, an unambiguous explanation for these trends is diffi¬ 
cult to reach, even in presence of slopes in the A[X/H]—Tc 


“ IDL (Interactive Data Language) is a registered trademark of 
Exelis Visual Information Solutions. 


relation steeper than Galactic chemical evolution effects, as 
they ma y not indicate necessarily the prese nce or absence of 
planets. iGonzalez Hernandez et ahl (l2013ll have found that 
after removing the Galactic chemical evolution effects in a 
sample of 61 single FG-type main-sequence objects, stars 
both with and without planets show similar mean abun¬ 
dance patterns. On the other hand. Galactic chemical evolu¬ 
tion is not capable of producing specific element-by-element 
difference between components of a binary system, whereas 
the plane tary accretion sce nario might reproduce it natu¬ 
rally (see IMack et al.ir2014ll . 

5. Origin of the elemental abundances difference 

The two stars have distinguishably different chemical com¬ 
position (see Fig. ISl), with a mean elemental abundance dif¬ 
ference between the N and the S components of ±0.067 ± 
0.032 dex, and statistically significant trends with the con¬ 
densation temperature (Fig.|4]). It is normally assumed that 
individual components in binary systems share the same 
origin and initial bulk abundances as they form within 
the same cloud. Therefore, under the reasonable assump¬ 
tion that the XO-2 components were born together from 
the same molecular cloud, their dissimilar elemental abun¬ 
dances seem to favor other explanations compared to Galac¬ 
tic chemical evolution and/or stellar birthplace effects as 
poss ible reasons of the abun dance patterns we observe (see, 
e.g.. lAdibekvan et al.ll20lll . 

Differences in elemental abundances between the XO-2 
system components could then reflect the likely different 
formation and evolution histories of the two planetary sys¬ 
tems orbiting around them, which might have led to metal 
deficiencies in XO-2S or excesses in XO-2N. We recall in 
fact the different architectures of the two planetary systems: 
XO-2N hosts a transiting hot Jupiter, with preliminary ev¬ 
idence for another object in an outer orbit (see Paper H, 
and references therein); XO-2S hosts two giant planetary 
companions on moderately close-in orbits (and possibly a 
third long-period object; see Paper I). 

We estimate the amounts of heavy elements being 
accreted or depleted in one of the c ompo nents follow- 
ing the prescriptions of [Murra^^e^j^yjl^OW L as done in 
iDesidera et al.l (|20n4ll and Desidera et al.1 (|2006ll . The mass 
of the convective envelope results to be about 0.06 Mq. 
Gonsidering the stellar [Fe/H] and A[Fe/H] between the 
components from Table [2j it results that the difference of 
iron within the convective zone is of ~ 10 Mq. Assuming to 
first order a similar composition of meteoritic material, this 
corresponds to 50 Mq of heavy elements. These amounts 
are larger than those considered in similar studies on plan¬ 
etary accretion (e.g., Pinsonneault et al.ir200lL iMack et al.l 
120141 , iTucci Maia et al.ll2014f) . due to the size of the convec¬ 
tive envelope and the super-solar metallicity of the XO-2 
components. We consider in the following sub-sections the 
two alternatives of heavy elements depletion in XO-2S and 
heavy elements excess in XO-2N. 


5.1. Heavy elements depletion in XO-2S? 

The existence of the abundance difference between the com¬ 
ponents of the XO-2 system can be interpreted i n the c on- 
text of the scenario devised by I Melendez et al.l (120091) as 
depletion of heavy elements in the stellar atmosphere of 
XO-2S, due to their being locked in planetary companions. 
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Table 3. Mean internal errors in abundance determinations due 
to given uncertainties (ATeff, A logy, and A^) in stellar param¬ 
eters for XO-2N and XO-2S. 


XO-2N 

Tefi = 5290 K 

log 5 = 4.43 


= 0.86 km/s 

A 

ATeff = ±18 K 

Aloggr = ±0.10 

<1 

= ±0.06 km/s 

[C/H| 

0.02 

0.05 


0.02 

[N/H] 

0.01 

0.05 


0.01 

[O/H] 

0.00 

0.05 


0.02 

[Na/H] 

0.013 

0.003 


0.010 

[Mg/H] 

0.003 

0.002 


0.005 

[Al/H] 

0.013 

0.013 


0.010 

[Si/H] 

0.005 

0.001 


0.003 

[S/H] 

0.010 

0.006 


0.006 

[Ca/H] 

0.016 

0.002 


0.015 

[Sc/H] 

0.002 

0.003 


0.005 

[Tii/H] 

0.022 

0.002 


0.026 

[Tiii/H] 

0.003 

0.002 


0.016 

[V/H] 

0.024 

0.001 


0.008 

[Cri/H] 

0.017 

0.002 


0.016 

[Crii/H] 

0.009 

0.003 


0.013 

[Mn/H] 

0.017 

0.001 


0.029 

[Fel/H] 

0.010 

0.001 


0.015 

[Fell/H] 

0.013 

0.004 


0.014 

[Co/H] 

0.008 

0.002 


0.003 

[Ni/H] 

0.003 

0.001 


0.018 

[Cu/H] 

0.006 

0.003 


0.017 

[Zn/H] 

0.005 

0.001 


0.019 

[Y/H] 

0.003 

0.003 


0.030 

[Zr/H] 

0.001 

0.003 


0.023 

[Ba/H] 

0.006 

0.001 


0.043 

[La/H] 

0.003 

0.006 


0.003 

[Nd/H] 

0.004 

0.004 


0.007 

[Eu/H] 

0.001 

0.006 


0.001 

XO-2S 

Tefi = 5325 K 

log g = 4.420 


= 0.93 km/s 

A 

ATeff = ±37 K A log 5 = ±0.094 A^ 

= ±0.03 km/s 

[C/H| 

0.04 

0.05 


0.01 

[N/H] 

0.02 

0.05 


0.01 

[O/H] 

0.00 

0.05 


0.01 

[Na/H] 

0.023 

0.020 


0.003 

[Mg/H] 

0.013 

0.015 


0.002 

[Al/H] 

0.022 

0.007 


0.003 

[Si/H] 

0.011 

0.001 


0.004 

[S/H] 

0.020 

0.020 


0.001 

[Ca/H] 

0.030 

0.026 


0.009 

[Sc/H] 

0.006 

0.034 


0.004 

[Tii/H] 

0.043 

0.009 


0.013 

[Tiii/H] 

0.008 

0.029 


0.010 

[V/H] 

0.046 

0.006 


0.003 

[Cri/H] 

0.034 

0.016 


0.015 

[Crii/H] 

0.020 

0.029 


0.005 

[Mn/H] 

0.035 

0.015 


0.015 

[Fel/H] 

0.020 

0.011 


0.008 

[Feii/H] 

0.026 

0.035 


0.006 

[Co/H] 

0.017 

0.017 


0.003 

[Ni/H] 

0.006 

0.004 


0.009 

[Cu/H] 

0.016 

0.006 


0.006 

[Zn/H] 

0.012 

0.005 


0.009 

[Y/H] 

0.002 

0.028 


0.015 

[Zr/H] 

0.002 

0.037 


0.009 

[Ba/H] 

0.006 

0.015 


0.020 

[La/H] 

0.005 

0.040 


0.001 

[Nd/H] 

0.005 

0.039 


0.003 

[Eu/H] 

0.001 

0.040 


0.001 


More precisely, our res ult appears qualitatively similar to 
that recently found by iTucci Mala et all (1201411 in the bi¬ 
nary system 16 CygA/B, in which the secondary compo¬ 
nent, which hosts a giant planet in eccentric orbit, is metal 
deficient with respect to the companion, which is not known 
to have planetary companions. Moreover, as found by the 
same authors, here we reveal a hint of the break in Tq 
between volatil es and refractories at ^ 800 K (Fig. d]). As 
pointed out by iTucci Maia et all (l2014fl , this means that 
the rocky core of the giant planet most probably was not 
formed in the inner disk, but at larger distances where they 
are more likely to form. 

The large amounts of heavy elements to be locked in 
planetary companions is challenging to estimate. However, 
it is a possibility that can not be dismissed in a super 
metal-rich environment such as that of the XO-2 system. 
Indeed, there are cases of compact transiting giant planets 
whose high density appears to be due to a very large con- 
tent of heavy elements, such as the ho t-Saturn HD149026b 
(jSato et al.H2005l iFortnev et al.ll2n(i^ . 

Theoretical models also predict the existence of such 
high-density Sat urn-mass planets form ed by up to 80% of 
heavy elements (iMordasini et al.ll2012f) . however they are 
most likely to form in th e outer region of a super metal-rich 
proto planetary disk (iDodson-Robinson fc Bodenheimerl 
[2onm . so that substantial migration would have been 
needed to carry such a planet at the observed separation 
of XO-2Sb. 

The inner planet XO-2Sb has a mass similar to 
HD 149026b and stellar metallicities are also comparable. 
Unfortunately, search for transits of XO-2Sb have been in¬ 
conclusive up to now (APACHE Team, priv. comm.), pre¬ 
venting the determination of the planet density. The pres¬ 
ence in XO-2S of two planets (with XO-2Sc more massive 
than XO-2Nb) and the possibility of an additional object 
responsible for the long term RV trend could further con¬ 
tribute to the segregation of heavy elements by the planets, 
favoring the locking scenario. The moderate eccentricities 
of the XO-2Sb and XO-2Sc (e ~ 0.2) could imply a smooth 
migration history of the known giant planets around XO- 
2S, possibly different from a dynamical instability event 
that might have caused the inward migration of XO-2Nb 
to its current tight orbit. Such a smooth evolution could 
also have allowed formation and survival of additional low- 
mass planets that remain undetected in the available RV 
data. This possibility was also mentioned by iTeske et al.l 

Iml). 

5.2. Heavy elements excess in XO-2N7 

The enrichment of refractory elements in XO-2N when com¬ 
pared with XO-2S could be due to accretion of material 
onto the star and subsequent ingestion of rocky material 
coming from already formed planets or from the inner part 
of the disk and pushed into XO-2N by its hot Jupiter as it 
migrated inward to the present p osition, most likely as a re¬ 
sult of planet-planet scatte ring (iGonzaled 119971 : llda &: LinI 
I20fl8l : iRavmond et al.ll20flh . Indeed in Paper H, it was dis¬ 
cussed that migration to the current close orbit of XO-2Nb 
can not be due to Kozai interaction for the present orbit 
configuration, supporting a phase of dynamical instability 
in the planetary system. The almost null eccentricity and 
the low spin-orbit angle of XO-2Nb (see Paper H, and ref¬ 
erences therein) are not relevant for inferences on the origi- 
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nal system configuration and migration history, as they are 
expected to be due to tidal effects considering the planet 
characteristics, the size of the convective envelope of the 
star, and the system age. 

The dynamical instability could have been triggered 
by modifications of the binary orbit caused, e.g., by stel¬ 
lar encounter in a dense stellar environment or a former 
hierarchical triple system that became unstable after the 
planets completed their g rowth a round the primar y star 
i|Marzari k, Barbierill2007() . Indeed. iKaib et all (1201311 stud¬ 
ied very wide binaries, indicating that they may reshape 
the planetary systems they host, changing the orbital ec¬ 
centricity or causing the ejection of planets or favoring the 
ingestion of material by the host star. In particular, they 
demonstrated that in binary systems with stellar masses 
and mean separation similar to those of XO-2N/S, 90% of 
instabilities occur after ^ 100 Myr, i.e. after planet forma¬ 
tion is complete. This might cause ingestion of planetesi- 
mals and/or planet ejections of one or more objects. This 
scenario has the advantage that the dynamical instability 
and the successive system reconfiguration could have oc¬ 
curred well after the shrinking of the outer convective zone 
of the star after the pre-main sequence phase, then avoiding 
one of the difficulties in the planetary accretion scenario. 

A challenge can be represented by the different situation 
of the HAT-P-1 system, whose secondary also hosts a hot- 
Jupit er with mass of 0. 53M.T and orbital period 4.46529 
days (|Bakos et al.l l2007ll . so not markedly different from 
XO-2Nb. The masses of the components of the system are 
slightly larger than those of the XO-2 system (I.I6M0 
and I.I2M0) and the binary separation is also large (1550 
AU) as for the XO-2 system. However, no significant abun¬ 
dance diff erence between t he components have been re¬ 
ported by iLiu et ahl (l2014f) . It was suggested that while 
the smaller size of the convective envelope could make eas¬ 
ier the imprinting of any planet signature in the system, 
the larger stellar mass favored a faster system evolution, 
with accretion before the shrinking of the convective zone, 
as more massive stars s eem to have shorter disk lifetime 
(|Williams fc Ciez'^l2011ll . More likely, the stochastic nature 
of the dynamical instability events should imply different 
outcomes for what concerns the ingestion of already formed 
planets and planetary debris. 

The amount of heavy elements to be ingested in the 
XO-2N enrichment scenario obviously corresponds to that 
assumed to be locked into planetary bodies in the XO-2S 
depletion scenario. As shown above, low-mass giant plan¬ 
ets could have heavy elements content of several tens of 
Earth masses, making the ingestion of an individual planet 
with peculiar characteristics a possible cause for the whole 
abundance difference. 


6. Conclusions 

We performed an accurate differential analysis of abun¬ 
dances of iron-peak, a-, s-process, r-process, and other ele¬ 
ments for the dual-planet-hosting binary XO-2, in which the 
stars show very similar stellar parameters (Tea ^ 5300 K, 
log5 ~ 4.43, M* ^ O.97M0, R^, ^ l.Oi?0, and [Fe/H]~ 
0.35 dex). To our knowledge, this is the first binary system 
in which both components host planets and whe re elemental 
abundance differences we re revealed (see also iTeske et al.l 
I2OI5I : iRamirez et al.l[20THl . 


We find significative abundance differences in almost all 
25 elements, in the sense that the slightly cooler compo¬ 
nent XO-2N is also the star with higher content of elemen¬ 
tal abundances, with stronger enhancement in the refracto¬ 
ries. The abundance differences range from ~ 0 dex for the 
volatiles up to ~ 0.11 dex for the refractory elements, with 
a significance at the 2a level for almost all elements. 

The differential abundances of XO-2N relative XO-2S 
show positive trends with Tq, with a significance of the 
slope at the 5 (t level considering all elements. This could 
mean that both components experienced the accretion of 
rocky planetary interior to the subsequent giant planets. 
We discuss about the possibilities that this result could be 
interpreted as the signature of abundances deficiencies in 
XO-2S or excesses in XO-2N, i.e. due to the ingestion of ma¬ 
terial by XO-2N or depletion by locking of heavy elements 
by the XO-2S planets, respectively. We estimate a mass of 
several tens of Earth as amount in heavy elements being ac¬ 
creted/depleted. For the time being, we do not know which 
of these two scenarios dominates or which one is the most 
appropriate and conclude that both are plausible. Future 
detailed numerical simulations and theoretical interpreta¬ 
tions will be needed to try to disentangle the skein and 
clarify whether one of the two hyphoteses is most proba¬ 
ble or if both are at work in such system. This is out of 
the aim of this work, which was finalized to look at possi¬ 
ble differences in elemental abundances in a planet hosting 
binary. 

We confirm the potentiality of the binary hosting stars 
as useful tools to derive in a very accurate way elemental 
abundance differences. Further detailed studies on chemical 
composition of double systems hosting planets can provide 
important constraints on stellar/planetary formation and 
subsequent evolution. Studies on how planet formation may 
influence the chemical composition of hosting stars can be 
important to select future exoplanet searches. 
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